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Abstract 

Contraction theory is a recently developed dynamic analysis and nonlinear control system 
design tool based on an exact differential analysis of convergence. This paper extends contrac- 
tion theory to local and global stability analysis of important classes of nonlinear distributed 
dynamics, such as convection-diffusion-reaction processes, Lagrangian and Hamilton-Jacobi 
dynamics, and optimal controllers and observers. 

The Hamilton- Jacobi-Bellman controller and a similar optimal nonlinear observer design 
are studied. Explicit stability conditions are given, which extend the well-known conditions 
on controllability and observability Grammians for linear time-varying systems. 

Stability of the Hamilton-Jacobi dynamics is assessed by evaluating the Hessian of the 
system state along system trajectories. In contrast to stability proofs based on energy dissipa- 
tion, this principle allows to conclude on stability of energy-based systems that are excited by 
time-varying inputs. In this context, contraction can be regarded as describing new variational 
conservation laws and the stability of entropy producing processes. 
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1 Introduction 



This paper shows how the tools of contraction theory can contribute to the stability and conver- 
gence analysis of nonlinear partial differential processes |I5J |71 [T^ [T51|. Contraction theory is a 
recently developed analysis and control system design tool based on an exact differential analysis 
of convergence. While differential approximation is the basis of all linearized stability analysis, 
what is new in contraction theory is that differential stability analysis can be made exact, and in 
turn yield global results on the nonlinear system ifTTlfTSl . 

Specifically, the paper considers stability and convergence analysis of vector nonlinear dif- 
fusion processes and Hamiltonian dynamics, and more broadly of nonlinear partial differential 
equations in the general form 

^ + h,{<l>, V0„ x,t) = V- G,,(V<I', X, t) (1) 
at 

where $ is the n-dimensional continuous state vector (whose components (pi consist, for instance, 
of chemical concentrations or mechanical velocities), and x is the m-dimensional vector of under- 
lying coordinates. Specific application examples are then discussed, including distributed chemical 
processes, classical fluid systems such as Bernoulli and Navier-Stokes dynamics, Lagrangian sys- 
tems, Hamilton-Jacobi dynamics, and Hamilton-Jacobi-Bellman designs of optimal controllers or 
observers IDSHHl. We also show that approximations of such processes using basis functions are 
themselves contracting, regardless of the coarseness of the approximation. 

Basically, a nonlinear dynamic system is called contracting if initial conditions or temporary 
disturbances are forgotten exponentially fast, i.e., if trajectories of the perturbed system return to 
their nominal behavior with an exponential convergence rate. In many physical systems this behav- 
ior can be interpreted as entropy producing — when the state corresponds to a physical quantity as 
velocity, temperature or a chemical concentration, the variation between actual and nominal behav- 
ior uses up available energy. As shown in LiSJ (to which the reader is referred for more details), it 
turns out that, for ordinary differential equations, relatively simple conditions can be given for this 
stability-like property to be verified, and furthermore that this property is preserved through basic 
system combinations, such as parallel combinations, series or hierarchies, and specific feedback 
combinations. 

In order to extend these results to the distributed case ([T]), we will regard any distributed quan- 
tity in the following as a limit of a spatial discretization. Accordingly the V-operator will be 
regarded as a limit of an explicit differentation and the / -operator as a limit of a sum. Particular 
instances of these extensions were obtained in lfT9tl20ll . and in lfT3lfT4ll using operator notation, 
which we adopt here. In the following V- corresponds to a inner product with the V operator, 
while V represents an outer product with the V operator. 

Section |21 analyzes the contraction behavior of nonlinear first-order partial differential equa- 
tions and their numerical solution, and discusses specific applications. Section |3] studies similarly 
second-order partial differential equations. Brief concluding remarks are offered in section |4l 
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2 First-order partial differential equations 



This section concentrates on first-order nonlinear partial differential equations on a continuum V, 
such as e.g. Hamilton- Jacobi equations, in the general form 

^ + /i,(<l>,V0,,f,t) =0 (2) 
ot 

Motivated by contraction analysis of ordinary systems in ifTSl . consider two neighboring so- 
lutions of Q at fixed time t, and the virtual displacement 5$. This leads to the virtual dynamics 

In most physical and technical applications, such as Bernoulli, Navier-Stokes and Hamilton- Jacobi 
systems, the term corresponds to an underlying velocity field of the quantity (pi in @, and so 
we will refer to this term as the flow "velocity". 

Existence and uniqueness questions in nonlinear partial differential equations are notoriously 
difficult |71|5l|T3|. Because the solution of ([T]) or Q needs not be differentiable with respect to x, 
V0 has to be carefully defined. We choose to define V^i causally along flow lines, i.e., so that it 
corresponds, in the sense of distributions, to the limit of an explicit backward differentiation against 
the corresponding velocity component (in the case that a component is zero, we can choose 
the direction of the differentiation arbitrarily without affecting Q). Accordingly we assume (pi to 
be given at all infiowing boundaries, i.e., parts of the continuum boundary dV where n < 0, 
with n the unit vector along the outward normal. Interpreting © as an ordinary infinitesimal 
discretization in space can then lead to a natural existence condition on the solution of ([T]). 

In this section, we derive two main results, one for a fixed continuum, and another for a con- 
tinuum moving and deforming according to a predefined velocity field in x space. Applications 
of the first result includes transport laws and Bernoulli dynamics, while applications of the second 
includes Hamilton- Jacobi dynamics and optimal controller/observer design. 



2.1 First-order p.d.e.'s on a fixed continuum 



Let us analyze the contraction behavior of © on a fixed region V . We first derive the main result 

dhj 

We start with the flux equation, for a continuum V of boundary dV, 



using partial integration, under the assumption that S(f)i and are continuously differentiable, 
before extending it to possibly discontinuous 6(f)i and -j^j- using a more technical proof. 



V ■ ( S(l),^6(p) dV= [ 6(P,-^ n 50, ddV 
oV(t)i J Jqv dV(j)i 



3 



where n is the unit vector along the outward normal of dV, and n the corresponding outward 
flux. Assuming that (pi is given at all inflowing boundaries, and using 



the time derivative of jy 54>i6(pi dV can be computed with partial integration as 

dh- V dh- 1 r dh- 



This yields the bound and main result 



2dtjy - Jy \ d<!> 2 dV^J ~ Jv 



where Ay denotes the largest eigenvalue of the symmetric part of — |^ + y 



dh 

a* ' 2 av* 



Assume now that Ay is uniformly strictly negative (i.e., 3 /3 > 0, V$, Vt > 0, Ay < — /? < 
0) then from Q any Jy dV converges exponentially to zero. As in contraction analysis 

for ordinary differential equations, this immediately implies, by path integration, that any finite 
difference between two arbitrary solutions converges exponentially to zero. Thus, as in stable 
linear time-invariant systems, the initial conditions are exponentially "forgotten." 

In the case that 6 (pi or is not continuously differentiable, inequality ^ can still be ob- 
tained using a more technical proof, as discussed in the Appendix. 

The argumentation can be extended to the following cases: 

• Consider again the dynamics Q with an upper bounded — 1| and a given initial condition 
$o. This implies that remains equal to zero Vt > 0, which implies uniqueness of the 
solution Vt > 0. 

• Apply now to system Q a set of linear constraints in $ or V$. Such a case may describe, 
for instance, mechanical systems with kinematic constraints, incompressible fluid flows, and 
chemical systems in partial equilibrium. The constrained dynamic equations take the form 

^ + /ii($,V0i,f, t)+pi = Q 

where the constraint terms p($, V$, x, t) are orthogonal to the constraint plane. Interpreting 
the above as an ordinary limit of a spatial discretization, we can conclude from 1.1 8J that 
contraction is preserved. 



In summary. 
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Theorem 1 Consider the dynamics 

^ + /i,(<l>,V0i,f, = 

with time t, m-dimensional coordinates x, n-dimensional state vector $, p constraint terms p 
orthogonal to p linear algebraic constraints in $ and V$, and V in line i defined as the limit of a 
backward differentation against the velocity 

For given (pi over all inflowing boundaries dV ( i. e. n < j and uniformly negative deflnite 
— H + I V • the system converges exponentially to a single solution, independent of initial 
conditions. The convergence rate is | Av|, where Ay is the largest eigenvalue of the symmetric part 
of-^ + i V ■ 



dh 



In the autonomous case ( hi = hi{^, V(pi, x), and constant boundary conditions and algebraic 
constraints ) and under the same conditions, the system converges exponentially to a unique steady - 



state ^ = 0. 



We will call a system contracting if Theorem [T] applies. Similarly, we call the system semi- 
contracting for -II + i V ■ < and indifferent for -|| + i V ■ = 0. 



More general local versions of the above theorem can also be derived as in 1 1 8 1 . Note that as the 
examples illustrate, contraction behavior can be interpreted physically as entropy producing, since 
5(i) could have been used to spend available energy instead, if the state corresponds to a physical 
quantity such as velocity, temperature, or chemical concentration. 

Example 2.1 : Consider the conservation of a quantity (j) in an m-dimensional continuum V 

d(h 

-^ + v{x,t) V0 = O 

with (/) given at all inflowing boundaries. Computing 

dh 1 ^ dh 1 ^ 
\ V • = -V -v 

8(1) 2 dV(t) 2 

we conclude with Theorem^on contraction behavior for a compressing flow field with V • w < 0, with 
minimal convergence rate | | . 

Consider another typical conservation of (j) in an m-dimensional continuum V 

^ + v-(Mx,t)) = o 



with (/) given at all inflowing boundaries and equilibrium points. Similarly to the above we can conclude 
on contraction behavior for a expanding flow field with V • u > and minimal convergence rate | | . 

Note that in the case that the velocity is itself the output of a contracting system, as e.g. will be the case 
later in example 13.31 then the overall system constitutes a hierarchy of contracting systems, and as such 
is contracting il8ll . 

Observer designs for conservation dynamics will be illustrated in example l331 and example IT6l D 
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Example 2.2 : Consider the m-dimensional Bernoulli dynamics 



^ + i V(l)^H-\x, t)V(t> + U{x, t)+p = 

with kinetic energy ^ Vcjr'" H~^{x, t)V(t>, potential energy U {x, t), and (j) given at all inflowing boundary 
conditions of the continuum. The constraint term p (energy due to pressure) stems from the incompress- 
ibihty condition 

Vh ■ H-^V(p = 

where V// denotes the covariant derivative with respect to the metric H{x,t). This dynamics governs 
many physical phenomena in fluids, structural mechanics and electromagnetism. Since 

_dh 1 - n 

the Bernoulli dynamics is indifferent. D 



Example 2.3 : Approximation of continuous distributed processes is routinely performed for numerical 
simulation, as well as in machine learning and vision. Let us now briefly discuss approximating the state 
of a contracting system Q as <l>(x, t) = ^{x, t)a, using n basis functions ^{x, t) and their coefficient 
vector a. The approximation leads to 

— + h($,V$,f,t) + p = (5) 

where the constraint terms p($, V$, x, t) refers to the error introduced by the approximation. Multiply- 
ing ^ from the left with and integrating results in an udpate law for a, 

j w^w dV k + y ^h(wa, Vw a, x,t) + ^ a^j dV = (6) 

where we used the orthogonality condition jy w^pdl/ = 0. The constraint terms can be computed by 
solving ^ analytically with t) = w(x, i)a and 

While the above is standard, the key remark is that, since a basis function approximation ^{x,t) = 
^{x, t)a can be interpreted as a linear constraint, contraction is preserved by the approximation (using 
the basic result in 1 18 1 on adding linear constraints to contracting systems). Note that the approximation 
technique can be extended to more general partial differential equations, with the same argumentation on 
the contraction behavior. 

Also note that we can smoothly extend the basis functions set with additional basis functions 'Wi{x, t) at 
any time during the simulation, by initializing the corresponding coefficient to = 0. Conversely, the 
deletion of a basis function Wi{x, t) from the set will lead to an additional disturbance Wi{x, t)ai and can 
hence only be accepted for small \ai\. D 

Finally, note that Theorem[T]can be extended to the following cases: 
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• Consider a spatial discretization of V. Then Theorem [T] can still be applied with a corre- 
sponding discretized V-operator since the method of proof in the Appendix is based on such 
a discretization. 

• Assume that the V-operator is generalized from a Jacobian to a covariant derivative with 
respect to a symmetric metric M(x, t). Within a local normal coordinate system fP\ defined 
by the metric, the definition of V in Q is unchanged (vanishing Christoffel term), so that 
Theorem [T] can be applied unchanged. 



2.2 First-order p.d.e.'s along flow lines 

Whereas the convergence of (p is important for p.d.e.'s on a fixed ontinuum the convergence of x is 
important when we consider the dynamics along flow lines e.g. for mechanical systems or optimal 
control. 

To assess the contraction behavior of x let us compute the evolution of the Hessian of along 
the flow lines. We will see later on in the optimal observer and controller dynamics and for physical 
Hamiltonian systems that the definiteness of VV0 will play an important role in the contraction 
behavior of the ordinary plant dynamics, using ifTSll . 

Let us compute the first and second gradient of © 

aV0 dh dh 



dt dx dVcj) 



dt dap dxdS/cf) dVcjydx dVcj)^ dWcj) 

Following system trajectories 

■ * (7) 



dV(t> 



leads to 



(iV(/) _ dVcf) ^dV(t)dx _ _dh 

dt dt dx dt dx 

dVVcf) _ dVV4> dVV(f)dx 

dt dt dx dt 

d^h _ d^h ^^^r d^h _ yy^T^!^yy^ 

dcP dxdV(j) dV(j)dx dVcjP' 

where the above is known as the ordinary momentum dynamics of a Hamiltonian (see e.g. [7J). 
Equation Q generalizes the Riccati equation (see e.g. ||4||) to nonlinear Hamiltonian systems. 

Let us perform a modal decomposition of Wcf) = AX with real eigenvalue matrix A and 
orthonormal eigenvector matrix X. Equation Q then takes the form 

= .Ar - X (-^X- + X- V - A f + x) X- - AX^X-A 

dt dx^ \dxdV^ J V dV(j)dx J dVcj)^ 
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where the diagonal corresponds to n decoupled scalar Riccati equations of the eigenvalues A in A. 

Uniform positive definiteness of A (i.e. VV0 or convexity of 0) can be shown Vt > with 
given initial uniformly positive definite Ao or V V0o if A moves away from in the neighborhood 
A = and always stays upper bounded. 

Using the modular Riccati equation the first condition (moves away from 0) is verified if the 
diagonal elements 

^d^agUu^^lA (10) 



dp A=o \ doP 

are larger or equal to for j = I and become uniformly positive in one of the higher derivatives 
(j > 1) if the above diagonal elements are not uniformly positive already. 

This is equivalent to requiring that is positve semi-definite and the remaining nuUspace 

becomes uniformly positive in one of the higher derivatives (j > 1) before it eventually becomes 
negative. The Lie-derivatives are here defined as 



dx^ dt \ dx^ J dxdV(f) dx^ dx^ dS/(j)dx 

The inverse Riccati dynamics is given as 

^Xl^ ^ VV0-^— VV0-^ + VV0-^— — + ^VV0-^ + (11) 

dt daP dxdV(j) dV(j)dx dS/<fP 

Performing similarly to the above a modular decomposition, the second condition (boundedness) 
is verified if the diagonal elements 

^^^^ = Aa,(x-L^g^xj (12) 

are larger or equal to for j = I and become uniformly positive in one of the higher derivatives 
(j > 1) if the above diagonal elements are not uniformly positive already. 

This is equivalent to requiring that is positive semi-definite and the remaining nuUspace 

becomes uniformly positive in one of the higher derivatives (j > 1) before it eventually becomes 
negative. The Lie-derivatives are here defined as 



_ d d^h \ ^ dh d^h ^ d^h dh 



The following theorem summarizes this conservation law of the curvature of along system 
trajectories. 
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Theorem 2 Consider the dynamics 



with time t, m- dimensional coordinates x, and action (p. Along trajectories, which move with flow 

dh 
9V0 



velocity x = -Mr , the Hessian W(j) evolves as 



dt dap dxdVcf) dV(t)dx 9V0^ 

An initial convex 4> (i-e. VV0o > Oj remains convex if conditions UUti and (17^ are fulfllled. 



Note that the conditions generalize the standard controllability Grammian condition ifTTl to general 
nonlinear Hamiltonian systems. Furthermore, they are explicitly computable, in contrast to the 
controllability Grammian which requires the unknown transition matrix. 



Example 2.4 : Consider the plant 

§ = /(x,n,t) 

with m-dimensional state x, j3-dimensional control input u, and cost-to-go dynamics 

-^ = -e{x,u,t) <0 (13) 

along the plant trajectory with final cost 4>f{x) at time tj. A typical example is £ = ^{x — x^)^ R{x — 
Xd) + \{u — Ud)^Q{u — Ud) , with desired trajectory Xd{t), supporting open-loop input Ud, and constant 
uniformly positive definite R and Q. (j) can be minimized locally by minimizing the Hamiltonian 

h = i{x, u, t) + f{x, u, t) 
over n leading to J| = V4> This leads to the well known Hamilton- Jacobi-Bellman equation l4ll 

^ + h{Vcj),x,t) = 

dx dh 



dt dVcf) 

on a continuum V which moves with This equation is solved backwards in time from final time 
tf. Assuming 4>f{tf) to be convex Theorem [fallows to compute under which condition stays convex 
Vt < t/. As a result Theorem Sallows to conclude when a global minium of (p is given by (fT3l . 

The ordinary contraction behavior of this optimal control design can be computed with the variation of 
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and using Q and STU 



— ((5x^VV0f) = -6x^ i —- W(P——tW<P 5x < 5x^L^—5x 
at ^ ' \ ox o\/(p^ ) ox 

We can conclude on (semi)-contraction behavior (i.e. unifomily positive definite W<j) and negative 
(semi)-definite L^^) if the conditions in and are fulfilled for j = 1 (j > 1). Note that 
asymptotic convergence can even be guaranteed for semi-definite L^^q^ since the differential length 
cannot get stuck due to the higher Lie derivatives in ( ITUl . 

In the case that the Hamiltonian dynamics is solved numerically, e.g. with basis functions or a neuronal 
network the above argumentation stays unchanged if the Hamiltonian is extended with the constraint term 
in©. 

If u cannot be computed analytically, but rather numerically through a convex minimization process, 
then the overall dynamics is still indifferent, since it then corresponds to a hierarchy lITSl composed of a 
contracting convex minimization process and an indifferent Hamiltonian dynamics. D 



Example 2.5 : Consider again the plant 

with m-dimensional state x, p-dimensional measurement feedback u, and cost dynamics 

^ =^(x,n,t) > 

with initial cost (poi^) at time to- A typical example is £ = ^(y — ym)^ R{y — Vm) + \ ifQu , with 
measurement ym, measurement estimate y{x, t), and constant uniformly positive definite R and Q. cj) can 
be minimized by minimizing a shghtly different Hamiltonian 

h = —i{x, u, t) + f{x, u, t) 
over u, leading to J| = Vcj) This leads again to the Hamilton-Jacobi equation 

^ + h{V<P,x,t)=0 (14) 

which has to be integrated now with t. This is a generalization of the well-known Hamilton-Jacobi con- 
troller design to nonlinear optimal observer design, with the optimal solution obtained at the minimum, 
where = 0. 

The convexity of (f) and hence the global uniqueness of the mimimum of (p can be shown with Theorem 
|2l This result also hold for an approximated system if the constraint term is added to the Hamiltonian. 

Les us now show the relation of this observer design to the well-known Kalman observer (see e.g. I.11J ) 
for linear time-varying systems. The optimal state estimate x can the be found at the minimum of (j), i.e. 
for V0 = 0, leading with (fT4l to 

(iV0 _ dV(t> ^ dVcpdx _ dh dh ^ ^y^dx _ ^ 

dt dt dx dt dx dV(j) dt 
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This yields the following observer design in x 

dx dh ,„„ ,s ^ dh 

with V0 = and the inverse co variance matrix (or information matrix) VV^ defined in Theorem |2l □ 



3 Second-order partial differential equations 

This section discusses second-order nonlinear partial differential equations, such as e.g. diffusion 
equations on a closed continuum V, of the form 

^ = V-G(V$,f,t) (15) 

where $ is the n-dimensional continuous state vector, x is the m-dimensional vector of underlying 
coordinates, and t is time. 

In most physical and technical applications such as diffusion and Lagrangian dynamics, the first 

V operator and the V operator in the argument correspond to a limit of a left and right differentation 
or vice versa. Accordingly either $ or the projection of V$ needs to be given on dV, so that both 

V operators can be defined there. 

Interpreting (1131) under the above conditions as an ordinary infinitesimal discretization in space 
leads to a natural existence conditions for the solution of O. 

Similar to section |2l consider now two neighboring solutions of © at fixed time t, and the 
virtual displacement 5$ between them, leading to the virtual dynamics 

d „ dG „^ 

— 5^ = V- — — (5V$ 
dt 9V$ 

The time-derivative of Jy dV can be computed as 

14- f 5«<l>dV = - f SV<^^^6V^dV 



2dt Jy Jy 9V$ 

where we have exploited that the transpose of a left-differentation corresponds to a negative right 
differentation. Lower bounding > J^iju > with kij^i = for ij ^ kl allows to conclude 

on semi-contraction behavior independent of the boundary conditions. 

For given $ on dV the minimal contraction rate can be computed with a Fourier expansion as 

11 



where is the maximum length of the continuum along the i axis. By path integration, this 
immediately implies that any finite difference between two arbitrary solutions has an equivalent 
contraction behavior. Combining this result with Theorem[T]using the superposition principle 1 18 1 
leads to 



Theorem 3 Consider the dynamics 

^ + /ii(<l>,V0i,f,t)+p. = V-G,(V$,f,t) 

with time t, m-dimensional coordinates x, n-dimensional state vector $, p constraint terms p 
orthogonal to p linear algebraic constraints in $ and V$, and V in line i defined as the limit of a 
backward differentation against the velocity 



Let U be the diameter ( maximum length ) of the continuum along the i axis and lower bound 

dG 

the positive semi-definite > ^ijki > Ofora (perhaps time-varying) Dirichlet condition (i.e., 
specified on dV) and define AijM = for the remaining (perhaps time-varying) Neumann 
condition (i.e., V$(t) specified along the normal ofdV). 

For given (pi over all infiowing boundaries dV (i.e. n < Oj and uniformly negative 

definite — |^ + ^ V ■ the solution converges unique and exponentially to a single solution. 



a* ' 2 " av* 



independent of the initial conditions. The convergence rate is \\v — Yl^=i Y2T=i '''/a ^ \> where 
Xv is the largest eigenvalue of the symmetric part of — ^ + | V ■ 

In the autonomous case ( h = h($, V$,x), G = G(V$,a;) and with constant boundary 
conditions and algebraic constraints) and under the same conditions, the system converges expo- 
nentially to a unique steady-state ^ = 0. 

We will call the system contracting if at least one of the two above conditions is fulfilled. We 
will call such systems semi-contracting if the above matrices are only semi-definite. And finally 
we will call them indifferent if the above matrices are skew-symmetric. 

The method of proof implies that all the results on contracting systems in ifTSll and exten- 
sions from section 121 can be extended to contracting reaction-convection-diffusion processes, with 
boundary conditions acting as additional inputs to the system. For instance, any autonomous con- 
tracting reaction-convection-diffusion process, when subjected to boundary conditions periodic in 
time, will tend exponentially to a periodic solution of the same period. The convergence is ro- 
bust to bounded or linearly increasing disturbances. Finally, the method of proof also implies that 
any regular spatial discretization of the above process is contracting as well, when the explicit 
differentaion is performed against the velocity. 

Again, the following examples will illustrate that contraction behavior can be intpreted physi- 
cally as entropy producing, since 5cj) could have been used to spend available energy instead. 



Example 3.1 : Consider a wafer disk, whose temperature is controlled with a continuous external light 
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source [Wj. The dynamic equations are 



^ + /i (<A^ - </.^) = V • g*{r) Vct> with r = ||V<A||. 

with > the wafer temperature, (po{t) the external temperature, h a radiation constant, and given (p at 
the boundary of the disk. Assume that g* >0 and ^^gj^ > then the Jacobian 

_ dg* V(t) Vc/) ^ d{g*r) Vcp ^ 

is positive semi-definite. A specific saturated diffusion term may be of the form g* = tanh(ar)/(ar) 
(with a a constant), for instance. 

Computing — + ^ V • = — 4 h cj)^ , we conclude with Theoreml^Jon contraction behavior with 
minimal convergence rate 4 h(p^. 

This result may be used to design observers for the system. Indeed, it means that an open-loop (identity) 
observer for the system will converges to the actual temperature distribution with minimal convergence 
rate 4 hcj)'^. Similarly, assume now that the actual temperature is measured at the boundary of the disk. 
Using this measurement as a boundary condition for the observer state increases the convergence rate 
according to Theorem|3j while at the same time preserving consistency with the actual plant (i.e., keeping 
the actual temperature as a particular solution of the observer equations with their boundary condition). 
Such designs are illustrated in fl^ for a linear g. □ 



Example 3.2 : Consider the viscosity solution Q of the Hamilton- Jacobi equation 

^'^ + h{V(t),x,t) = g V^(/> 



dt 

aV0 



along a continuum V, which moves with the velocity -Mh- and given tp or V(/> along the normal of dV. 



Computing 

_dh 1 _dh_ _ 

we can conclude with Theorem |3l on semi-contraction for 5 < 0. The minimal contraction rate of the 
dynamics is governed solely by the viscosity term, the geometry and the boundary conditions of the 
continuum V . d 



Example 3.3 : Consider the m-dimensional Navier-Stokes equation 

+ Vvv + Vf/(x, t) + Vp = 5 V • Vv 

with velocity v, potential energy U{x, t), viscosity g, incompressiblity condition V • v = leading to a 
pressure gradient Vp and v given at all boundaries. Computing 

9h 1 _ 5h 1 

\ V • = - Vv 

5v 2 5Vv 2 
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we can conclude with Theorem|51on contraction behavior for positive definite Vv + YliLi ^ where 
li is the diameter of the continuum in direction i. Accordingly, uniqueness of Navier-Stokes solutions can 
be concluded for bounded Vv. 

Comparing this result with Euler dynamics (i.e. 5 = 0) we see that the stability is increased by the vis- 
cosity term. Comparing this result with the velocity potential (j) in the Bernoulli dynamics in example l2!2l 
potential instabilities of the Navier-Stokes can only be explained with vortices within the dynamics. This 
nonlinear result on a finite convergence region generalizes the well-known stability results on continu- 
ously differentiable infinitesimal small disturbances in 0. □ 



Example 3.4: Consider a continuum V with $ or given along the normal of dV, and the Lagrangian 

^ = T{V^,x,t) - U{<^,x,t) 

Assume that ^ is convex, i.e. that |^ and — are both uniformly positive definite. Then can 

be minimized using the Euler-Lagrange dynamics (see e.g. |4J) 



dt dV^ 
Computing 

_dh_ 1 _dh_ _ dHl_ 
contraction can be concluded using Theorem |3l 

This result is especially useful in function approximation and image processing, where — U describes 
e.g. the quadratic deviation between measured and estimated and T describes a nonlinear smoothing 
cost, e.g. penalizing small and medium gradients but letting through large gradients such as edges. □ 



Example 3.5 : In machine vision, consider an image (grey level) cj), with cf) given at all inflowing bound- 
aries of the imaging plane, and given along the normal of the remaining boundaries. Consider the 
following combination of optical flow (conservation of brightness) from example 12.11 and Lagrangian 
dynamics from example f3.4l (see 1 ,22. .23J ). 

with given camera motion flow v{x, t). Computing 

dh I dh _ d"^^ V • V 

we can conclude with Theorem |3lon contraction behavior for positive definite and strictly negative 

"3^ 2 • '-' 
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The approach can be vastly extended by allowing for a prior differential coordinate transforma- 
tion in $, as in |18 1. Specifically, the line vector 5^ between two neighboring trajectories can also 
be expressed using the differential coordinate transformation 6^ = 6^ where 0($, V$, x, t) 
is a square matrix. This leads to a generalization of our earlier definition of squared length 

/ 5^^d^dV= [ 6^'^M6^dV (16) 
Jv Jv 

with metric M = 0^0. Note that in general we cannot expect to find explicit new coordinates 
^{^,t), but 5^ and 6^^6^ can always be defined. We require M to be uniformly positive 
definite, so that exponential convergence of 6^ to also implies exponential convergence of 6^ 
to 0. The following example illustrates the application of Theorem |3] to new coordinates 5^. 

Example 3.6: Theorem |3]can be applied to a distributed version of an example studied in detail in 1 19 1. 
Specifically, consider the temperature-dependent reaction A —>■ B man open reaction volume V of size 
100 X 100 units 

with < Cyi < 1 the partial concentration of yl, T > the measured temperature, g > the diffusion 
constant, E > the specific activation energy, and v{x, t) the velocity field within the reaction chamber. 

The reaction volume is open at the left and the right and closed at the top and the bottom. Accordingly we 
assume that V = and V T = everywhere along the normal of the boundary, except for x = 
where we have Cyi = and T = 500. For 40 < y < 60 the chemical A is injected with concentration 
CAf{t) and temperature Tf{t). Extending |19L the observer can be defined similarly as 

i(?)--(^vi)--("v?)-(-:;„)---(::)/;-- 

An appropriate linear coordinate transformation with constant uniformly positive metric is computed in 
lIT^ for the ordinary reaction terms, based on a linear matrix inequality. Since the convection and diffu- 
sion terms do not change under such a linear coordinate transformation. Theorem |3] implies contraction 
behavior for both observer designs and the plant for > 0, V • w = 0, and uniformly negative definite 
projection of the Jacobian |^ on the new coordinates. 

The corresponding system responses can be accessed in |25 1. D 



4 Concluding Remarks 

This paper focused on the contraction analysis of important physical and engineered first and 
second-order distributed systems, such as Navier-Stokes, Euler, Lagrangian and Hamilton- Jacobi 
dynamics. In principle the technique can be extended to higher-order systems by repeating the 
partial integration. 
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Appendix: Derivation of inequality dJ) 



If either 5 (hi or 



dhj 



is not continuously differentiable, then standard partial integration cannot be 



used in the derivation of inequality Q. The inequality can still be obtained using a more technical 
proof. Specifically, we compute the convection term in 



av, 



d ^ , dhi -^-^ 



dhi 



9$ 



av 



based on the limit of a spatial discretization, where V is defined against the velocity component 



dhj 



and lower bound this discretization matrix with the diagonal matrix -^f 



dhi 



, SO that we 



get inequality Q again. This lower bound is obtained using the following computation 



E 



1 dhi 1 


1 dhi 1 


I dhi 1 \ 




1 avj</.i^l 


idS7j(f>i^\- J 



I Q^I I 






/ 1| dhi I 

/ 2\dV■i4>^k^ + 

I dhi I 



lav 







> 





I dhj I 





1 I dhj I 








I dhi I 



a/ii I 



k,k-l..k+l 
\ 



J 



k-l..k+l,k-l...k+l 



I av 



a/i; 



2 lav 



- J 



k-l..k+l,k-l..k+l 



with discretization index k. 



dh. 



av 



3<Pik' 



max 



ahj 



,0), 



av 



max 



ahj 



,0), and where 
avj0ik' ^' 



^avj<^ik 

only the non-zero elements of the discretization matrices are shown. In the case that </)j is given 
either at A; + 1 or /c — 1, then the corresponding row and column must be deleted, without any 
impact on the definiteness. 
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